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Zitterbewegung of neutral relativistic particles propagating along a constant magnetic and/or 
electric field is studied. It is shown that spin Zitterbewegung, when superimposed on the Larmor 
precession frequency, leads to a beating pattern. The existence of a forbidden frequency of spin 
precession is predicted. Modifications of position and velocity Zitterbewegung due to lifted spin 
degeneracy manifested in the appearance of longitudinal and transversal Zitterbewegung, each with 
two Zitterbewegung frequencies and resulting beating patterns, are reported. 
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I. INTRODUCTION 

The Dirac equation unified special relativity with 
quantum mechanics for fermions. It was a step further 
beyond the Schrodinger equation due to its invariancc 
with respect to Lorentz boosts. Several early experimen- 
tal successes of it were the accurate prediction of the in- 
trinsic magnetic moment of the electron, the prediction 
of antiparticles and the correct description for the fine- 
structure of the spectrum of the hydrogen atom. These 
were based on a form of the Dirac equation which de- 
scribes relativistic electrons. There exists a similar form 
of it , which offers a phenomcnological description of some 
compound particles like the neutron and the proton in 
the relativistic domain. 

Despite its successes regarding the relativistic behav- 
ior of the electron, certain solutions of the Dirac equa- 
tion predict some unusual effects, which have not been 
experimentally measured so far for neither the electron, 
the neutron or any other particle [l[ . The best known of 
them are (i) Klein's paradox [l[ — the tunneling of a rela- 
tivistic electron into an infinite potential barrier, and (ii) 
Zitterbewegung (ZB) [l[ — helicoidal motion of the ex- 
pectation value of the position f and oscillatory behavior 
of the expectation values of the momentum p, the spin 
S and the orbital angular momentum of a free electron. 
Recently these two effects have gathered renewed atten- 
tion by the scientific community, thanks to several pro- 
posals for their experimental emulation in physical sys- 
tems like solid-state quantum wells [|[ , trapped ions 0- 
@ , graphene Q and optical waveguides [|[ . 

Most of the recent work on ZB is focused on electrons. 
In contrast here we consider the Dirac equation for neu- 
tral particles propagating relativistically in external mag- 
netic and/or electric fields [l| along the direction of the 
applied field. The Dirac equation [lj for these particles 
deviates from the Dirac equation for the electron in three 
aspects: (i) the term coupling the electrostatic potential 
to the charge of the particle dissapears; (ii) the terms 
coupling the vector potential to the charge and momen- 
tum are not present; (iii) a term coupling the anomalous 
magnetic dipole moment (MDM) to the magnetic field 
and a term coupling the eventual electric dipole moment 
(EDM) to the electrostatic field appear. 



II. ZB OF FREE PARTICLES 



The Dirac Hamiltonian for a neutral particle in the 
absence of external fields reads 



H = ca x -p x + (3m c 2 



(1) 



where c is the speed of light, uiq is the rest mass of the 
particle, a x and /3 are the Dirac matrices. The spectrum 
of the Hamiltonian is doubly degenerate and its eigenval- 
ues are E± = ±E, with E = \J (cp x ) 2 + (m c 2 ) 2 . It is 
well known that if the initial state is a superposition 
of both positive and negative energy eigenstates, in the 
process of evolution the expectation values of velocity, 
position and spin will oscillate with the same frequency 



,zb 



2E 



(2) 



This is the free-particle ZB. 

There exists a lower bound lj = 2mc 2 /U for the free- 
particle ZB of the three quantities S, a and r. In the 
literature [l(| it is usually accepted that the frequency of 
the free-particle ZB is approximately w = 2mc 2 /%. How- 
ever, for relativistic speed of propagation, cp x becomes 
comparable to mc 2 . Then there will be considerable con- 
tribution of cp x to the ZB frequency w zb , Eq. @. The 
difference 



+ (toqC 2 ) 2 



2mc 2 



(3) 



between the rest frame (p x = 0) and the observer frame 
of reference amounts to a blue shift of the ZB fre- 
quency: it is demonstrated in Fig. [TJ Therefore, the 
free-particle ZB frequency has the opposite behavior of 
the earlier reported [ll[ relativistic red shift effect for 
the frequency of Larmor precession. It is a consequence 
of the relativistic energy-momentum-mass relationship 
E 2 = (pc) 2 + (mc 2 ) 2 , and the fact that ZB is observed due 
to an interference between positive and negative energy 
solutions. 
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in the matrix elements between the eigenstates, 



0.2 0.4 0.6 0.5 
Velocity v (in units c) 



FIG. 1. Demonstration of blue shift effect of free particle 
ZB frequency uf h , Eq. ©. The ZB frequency ui zh increases 
with the particle velocity v, and respectively, its momentum 
p — jmv, 7 = 1/^/1 — (v/c) 2 . 
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Using the plane wave eigenstates \l,s) of the Hamilato- 
nian (j4j, an arbitrary wavepacket can be represented as 
the superposition 



jr( P -X-E?) 



(7) 



which includes both positive and negative energy eigen- 
vectors @. 



A. Modification to spin-ZB 



III. MODIFICATION OF ZB DUE TO LIFTED 
SPIN DEGENERACY 

Both Larmor precession of spin along an external 
field and spin-ZB manifest themselves as spin precession. 
However, despite their similarity there are important dif- 
ferencies between them. Spin-ZB is a consequence of in- 
terference between positive and negative energy compo- 
nents of the Dirac equation, while Larmor precession is 
a consequence of lifting of spin-degeneracy by external 
fields. It is natural then to expect some modifications of 
at least spin-ZB, and eventually position-ZB, when spin 
degeneracy is lifted by external fields in the Dirac equa- 
tion. 

We address these problems by examining a relativistic 
model of a neutral particle possessing MDM, and even- 
tually EDM, placed in an external field, in which the 
static magnetic and/or electric fields are applied along 
the propagation direction. Working with the a;-direction 
in SI units the Hamiltonian takes the form 

H = ca xPx + Pmc 2 + 2pS x (dE - fiB) , (4) 

where S x is the x-component of the spin vector operator 
in relativistic theory, E is the electric field, B is the mag- 
netic field, d is the size of the EDM and fi is the size of 
the MDM of the relativistic particle. The four distinct 
eigenvalues of the Hamiltonian are 



El = ±y / {cp) 2 + (mc 2 + A) 2 , (5a) 
Ei = ±^/(cp) 2 + (mc 2 - A) 2 , (5b) 

where A = dE — fiB. The four distinct eigenvalues of 
the Hamiltonian correspond to four different frequencies 



The expectation value (^IS^I^) of the component of 
the spin vector operator along the propagation direction 
does not contain an oscillating component. Its value 

(*|&|*) - ^(E ;=± E. 5=T ,Jc p , M | 2 ( S ,/|fe| S ,0) is a 

constant of motion, determined by the initial conditions. 
This is natural since in this geometry S x is the helicity 
of the particle. The expectation value of the y and z 
components of the spin operator read (j = y,z) 

0i) = ^E 2Re ( c ^,t c ^4(t^l^l^l)e^ t ) + 
7 P i=± 

+ fY, m << + ^-A^+\S ] \- 1 s)e^ t ). (8) 

They contain two oscillating terms with different frequen- 
cies: a relativistic Larmor precession frequency, 

ujl = ^ (cp) 2 + (mc 2 + A) 2 

-^M' + K-A) 2 , (9) 

and a modified spin-ZB frequency, 



ujf = -V(cp) 2 + (mc 2 + A) 2 

+ ^(cp) 2 + (mc 2 -A) 2 . (10) 

The presence of two different frequencies for every plane 
wave component leads to a beating pattern with an an- 
gular frequency 

LO Sh = uf -0J L = |VM 2 + (™C 2 - A) 2 . (11) 
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FIG. 2. Demonstration of blue shift effect for the spin-ZB fre- 
quency uj2 h , red shift for Larmor precession frequency ujl and 
blue shift effect for their difference - the spin beating fre- 
quency to Bh . The figure is plotted assuming A = 0.4 mc 2 



Earlier we have reported the existence of an up- 
per bound 2mc 2 /h for the Larmor precession frequency 
of a relativistic neutral particle. On the other hand, the 
lower bound for the free-particle spin-Zitterbwegung is 
not modified by the coupling of EDM and MDM to ex- 
ternal fields and therefore it is still 2mc 2 /h. However, in 
a realistic situation the two limits can not be achieved by 
lot, and w| b due to the Heisenberg uncertainty principle. 
Thus the frequency 



2mc 



^forbidden 



(12) 



is a forbidden frequency for spin-precession within the 
considered model. Below it, spin precession is due to 
Larmor precession, whereas above it, it is because of ZB. 
For the Hamiltonian ^ the transverse spin components 
precess with the two frequencies at the same time leading 
to the beating pattern with beating frequency uj sh . 

The modified ZB frequency w| b exhibits a relativis- 
tic blue-shift effect similar to the one discussed for free- 
ZB frequency and contrary to the Larmor precession 
frequency [llj, which is subjected to red-shift effect. The 
beating frequency uj sh also exhibits blue shift effect — it 
is larger for particles with larger velocity and respectively 
momentum. The behavior of the three frequencies versus 
the velocity is illustrated in Fig. [2J 



B. Modifications to velocity and position ZB 

While modification of spin-Zitterbwegung by MDM-B 
or EDM-E coupling is expected, the modification of or- 
bital dynamics by the coupling of an external magncto- 
static/clcctrostatic field to the spin-degree of freedom is 
more surprising. We investigate these by computing the 



expectation values of ctj and fj (j — x, y, z) with respect 
to the superposition ([7]). We first calculate the modified 
velocity-ZB in the Schrddinger representation and then 
the position-ZB from the expression 



|/*+<M)«(M)<i=* = c/*- 



(x,t)a x ^(x,t)d 3 x 
(13) 

and integration [9( of the expectation value of velocity 
with respect to time and multiplication by c. 

It turns out that Eq. (j4j leads to two different types 
of velocity and position ZB: (a) longitudinal for the x 
component of the velocity and position operators here, 
and (b) transverse for the y- and z- components of the 
velocity and position operators. They are characterized 
by different precession frequencies. 

The expectation value (a x ) for the longitudinal 
velocity-ZB reads 



(«x> 



E 

P Ls 



CPx 

Ei 



P ,i 

N 2 { Px )2Re 



c +,p,i c -,p,l e 



c +,p,t c ->P't e 



(14) 



V'<£t(£|+£j)(£t+-Ej) 



El = mc 2 



A. 



where Ni{p x ) 

and N 2 (p x ) = — r - 

The first term of Eq. (fLtj) describes the group veloc- 
ity of the wavepacket, obtained as the weighted sum 
of the averaged velocities of the components of the 
wavepacket (0. The other two terms describe the 
velocity-ZB for the Hamiltonian (QJ. Excluding the ini- 
tial position, the expectation value of the transverse 
velocity-ZB reads 



J P Ls P,l 



+ TN 2 {p x )ZRe 
7p 



7Vi( P:c )2Rc 



p 



: +,p4 c - p4 iu? t 



(15) 



where again the first term describes classical motion 
while the other two describe position-ZB. 

In contrast to the free particle case [l], Q where there 
is just one frequency for velocity- and position-ZB, in 
Eq. and Eq. (fT5"|) there arc two different frequencies 
of longitudinal velocity-ZB wf b = 2^/{cp) 2 + (mc 2 + A) 2 
and w§ b = 2y / (cp) 2 + (mc 2 — A) 2 . Both of them are 
di fferent from the free particle ZB frequency w zb = 
2\J (cp) 2 + (mc 2 ) 2 , by a term determined by the size of 
the spin-splitting A. In the rest frame of the particle 
they become 



Lof(0) = 

W f(0) 



2mc 

h 

2mc 2 



2A 

X 

2A 



(16a) 
(16b) 
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FIG. 3. Demonstration of blue shift effects fir longitudinal 
velocity and position ZB frequencies uif 3 and a>f b , and red 
shift effect for their difference - the orbital beating frequency 
uj° h . Plotted assuming A = 0.4 mc 2 . Clearly the longitudinal 
ZB frequency cuf" can take values below 2mc 2 /S, as well as 
the ZB beating frequency cj° b . 



Obviously, u>l is always greater than the lower limit- 
ing frequency of free-particle ZB 2mc 2 /h, while w| b can 
take values which arc smaller; this is illustrated in Fig. EI 
In principle, by application of strong enough magnetic 
field the longitudinal vclocity/position-ZB frequency can 
be lowered to the observable range for current technol- 
ogy. Thus when the spin degeneracy of the particle is 
lifted by an external field, 2mc 2 /ft is not a lower limit 
for velocity/position-ZB, while it is still a lower limit for 
spin-ZB, as shown above. 

Furthermore, the presence of two different frequen- 
cies in the expressions for longitudinal velocity/position- 
ZB leads to the appearance of a beating pattern with an 
orbital beating frequency 



wf = uf - u)f = 2^{cp) 2 + {mc 2 + A)< 
-2 y / {cp) 2 + (mc 2 -A) 2 



(17) 



The beating frequency is just twice the Larmor preces- 
sion frequency of spin-splitting u>° h = 2wl, which makes 
the predicted effect readily observable with current tech- 
nology. With respect to Lorentz transformations of the 
coordinate systems the two longitudinal velocity /position 
ZB frequencies exhibit blue shift effect similar to the one 
for the spin-ZB frequency. In contrast, the beating fre- 
quency, as their difference, manifests rclativistic red shift 
effect. This is the opposite behavior of the beating fre- 
quency w sb of spin-ZB. These dependences of the ZB fre- 
quencies on velocity are plotted on Fig. [3] 

The expectation values of the y and z components of 
velocity/position have similar structures as the expecta- 
tion values of spin-ZB. They have terms with two fre- 
quencies lol and w| b , which are different from the free- 
particle ZB and from the longitudinal velocity/position- 



ZB frequencies Lof 3 and ujf 3 . The position expectation 
values reads 

ft) = f p g 2Re[%^( t ,/|d J |U) e ^ 

+ i^2Re\ C ^f'2'-' S (s',+\a j \-,s)e 
J P s ^ s , L 1^2 



(18) 



The second term in Eq. (|18p . with the angular frequency 
w| b , is easily identified as position-ZB. It is the result of 
an overlap between the components of positive and nega- 
tive energies and of opposite spin/helicity clearly visible 
in the expectation value of ctj, j = y, z. The transverse 
position-ZB frequency wf b has the same properties, as 
the one discussed in the context of spin-ZB and, in par- 
ticular, it is subjected to the blue-shift effect. 

The first term in Eq. (fT8"|) with the angular frequency 
equal to the Larmor frequency ujj, is harder to be identi- 
fied as traditional ZB. While it describes a phenomenon 
similar to traditional ZB — oscillation in the expectation 
values of velocity and position, it does not stem from in- 
terference between solutions with positive and negative 
energy. As it is clear from the structure of the matrix 
element (f, i|dj|i, \) (I = ±) of the velocity operator ctj 
(j = y,z), it stems from matrix elements with the same 
sign of energy but with opposite spin orientations. Ex- 
plicitly, the matrix elements (t, i|aj|Z,4-) are given by 



(-, ;KI- t) + (+,t KI+4) = 

cp x (huj L + 2A) 



+ 



"/ 



_ cp x (fiuj L + 2A) 
V 



cp x (huj L - 2A) 

K 

(19a) 
cp x (huj L - 2 A) 



C 



(19b) 



where C = 2^j E^E 1 _{e1 + El){E l _ + E^) and rj = 

2y^|(4+4)C4+4)- 

Since £ ^ rj the amplitude of the oscillations will 
be dominated by the matrix elements (—,4- \oij\— A) 
(j = y,z) between the solutions with negative energy. 
Furthermore, the matrix elements depend on the size of 
the relativistic Larmor frequency wl , the interaction en- 
ergy A and the momentum of the particle. In the rest 
frame of the particle (p x = 0) the oscillations disappear. 
They also disappear when the interaction energy A be- 
tween the particle MDM/EDM and external static field 
is zero. This is a feature of relativistic spin splitting. The 
lifting of spin degeneracy leads to precessionary motion 
in the plane transverse to the direction of spin splitting. 
This shows that by lifting the spin degeneracy and by 
the orbital effect of spin splitting it might be possible to 
observe ZB with a frequency smaller than the limiting 
frequency 2mc 2 /% of a free particle. 
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The combination of two frequencies in the expectation 
values of (fj) leads to the appearance of a beating pattern 
with a beating frequency uj^ — w| fc — u>l- It has the same 
properties as the beating frequency of the modified spin- 
ZB. 

The lowering of the ZB frequency w| b below the free 
particle limit 2mc 2 /%, as well as the beating frequencies, 
lay down a route for eventual detection of ZB phenom- 
ena in conventional experiments. This is further helped 
by the red-shift effect for the longitudinal beating fre- 
quency w ob . An alternative route for experimental stud- 
ies of the above effects is quantum simulations. In a pre- 
vious work [l2j . we have proposed the emulation of the 
Dirac equation with EDM in trapped ions. The same 
scheme of simulation can be applied to the emulation of 
the effects considered here. After the preparation of an 
appropriate initial state the dynamics can be simulated 
in the proposed way [12]. The different ZB frequencies 
can then be measured using the technique of mapping 
spatial variables onto the internal states and subsequent 
measurements of them I6f . 



We predicted the existence of a "forbidden frequency" of 
spin precession in ID and showed that it is not modified 
by the presence of external fields. We also predicted a 
modification of the frequency of spin-ZB. We showed that 
the beating frequency stemming from the combination of 
spin-ZB and Larmor precession exhibits a motional blue- 
shift effect. Furthermore, we predicted the existence of 
two different position ZB phenomena: transverse and lon- 
gitudinal. The longitudinal position-ZB is characterized 
by two frequencies different from the free particle one 
2mc 2 /h and leading to a beating pattern. We showed 
that the two ZB frequencies exhibit a motional blue-shift 
effect, while the beating frequency exhibits a motional 
red-shift effect. The transverse ZB is also characterized 
by two frequencies and a beating pattern similar to the 
spin-ZB. However, while the higher frequency can be in- 
terpreted as ZB frequency, the lower can not be identi- 
fied as conventional ZB, because it stems from matrix 
elements with the same sign of energy but opposite spin 
orientations. We referred to it as orbital effect of spin 
splitting. 



IV. CONCLUSIONS 

In this work we pointed out the existence of a motional 
blue-shift effect for the ZB frequency of a free particle. 
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